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Weyl semimetals (WSM) are topologically protected three dimensional materials whose low energy excita¬ 
tions are linearly dispersing massless Dirac fermions, possessing a non-trivial Berry curvature. Using semi- 
classical Boltzmann dynamics in the relaxation time approximation for a lattice model of time reversal (TR) 
symmetry broken WSM, we compute both magnetic field dependent and anomalous contributions to the Nernst 
coefficient. In addition to the magnetic field dependent Nernst response, which is present in both Dirac and Weyl 
semimetals, we show that, contrary to previous reports, the TR-broken WSM also has an anomalous Nernst re¬ 
sponse due to a non-vanishing Berry curvature. We also compute the thermal conductivities of a WSM in the 
Nernst (VT _L B) and the longitudinal (VT || B) set-up and confirm from our lattice model that in the parallel 
set-up, the Wiedemann-Franz law is violated between the longitudinal thermal and electrical conductivities due 
to chiral anomaly. 


I. INTRODUCTION 


After the theoretical prediction of topological insulators, 
and their subsequent experimental realization, the field of 
topological condensed matter has grown manifolcP®. The 
topological order manifested in these systems is not associ¬ 
ated with spontaneous breaking of a symmetry, but rather can 
be described by a topological invariant which is insensitive to 
a smooth deformation of the Hamiltonian. Usually the robust 
topological protection is associated with a non-zero spectral 
gap in the bulk of the system, and the presence of protected 
zero energy surface states is regarded as the hallmark of a 
non-trivial topological phase of matter. However, recently it 
has been proposed that systems in three spatial dimensions, 
in the presence of broken time-reversal (TR) and/or space- 
inversion (SI) symmetry, can also be topologically protected 
even without a bulk energy gapP^l^. These are Weyl semimet¬ 
als (WSM) - the nomenclature based on the Dirac/Weyl equa- 
tion which is used to describe their low energy excitationP^l 

A number of recent experiments have claimed to be able 
to observe the Weyl semimetal phase in an inversion asym¬ 
metric compound TuAJ^^EII and also in a 3D double gyroid 
photonic crystal^, without breaking TR. Another route which 
can result in the experimental verification of the novel WSM 
phase is to first realize a 3D Dirac semimetal and then break 
time reversal symmetry by applying a magnetic field, which 
will split a Dirac cone into a pair of Weyl nodes. Na^B i and 
Cd 3 As 2 were recently proposed to be Dirac semimetalJ^HlII 
and also have been confirmed experimentally by a series of 
experiment^2211221 in Bii_a;Sb 2 ; for a; ~ 3 — 4% also the 
Dirac semimetal phase has been predictecP^®^, and experi¬ 
mental signatures of realizing a WSM phase by breaking TR 
have been reported^. 

A simple WSM with broken time reversal symmetry can 
be desribed by a pair of linearly dispersing massless Dirac 
fermions governed by the Hamiltonian: iT±(k) = fzhvpcr ■ 
(k — K±), where cr is the vector of Pauli spin matrices 
defined in the space of two non-degenerate energy bands, 
uf is the Fermi velocity, and K± are the two band touch¬ 
ing points separated from each other in momentum space by 
ko = K+ — K_. It is essential that kg is non-zero to ensure 


that the system breaks TR and is topologically non-trivial, in 
which case and H_ describe two Weyl fermions of op¬ 
posite chirality. The two band touching Weyl points act as a 
source and a sink (monopole and anti-monopole) of Berry cur¬ 
vature, which acts as a fictitious magnetic field on the electron 
wave-function in momentum spaced. For kg = 0, the two 
Weyl points collapse onto each other giving rise to a topo¬ 
logically trivial (i.e with a zero Berry curvature flux) mass¬ 
less degenerate Dirac fermion. The topological nature of a 
WSM leads to a host of interesting physics, for example Berry 
curvature induced anomalous transport, namely charge and 
ther mal H all conductivitie^^^Ms] and open Fermi arcs on sur- 
face^nEllEl Anomalous transport phenomena, however, have 
been already known to exist in a variety of systems whic h pos- 
sess a non-trivial distribution of the Berry curvature flutP^Ell 
In a WSM, more interestingly, each Weyl node is chiral, with 
the chirality quantum number protected by a quantized flux of 
the Berry curvature, also known as Chern flux, which results 
in another peculiar phenomen on known as chiral anomaly 
(or Adler-Bell-Jackiw anomalyj®^l2lllIl xhe chiral anomaly 
concerns with the nonconservation of chiral charge i.e. an 
imbalance of charge between two distinct species of chiral 
fermions in the presence of non-orthogonal applied electric 
and magnetic fields. Several transport signatures have been 
proposed to test chiral ano maly such as negative lo ngitudin al 
magenoresistancJi^I^IlllIIland chiral magnetic effecP^^^^E^of 
which the f ormer has been recently claimed to be observed in 
experimentJ^®®^. 


A Dirac node can be split into two Weyl nodes by breaking 
either the TR symmetry or SI symmetry. Figure [T] shows a 
linearly dispersing Dirac node split into a pair of Weyl nodes 
when TR symmetry is broken, and also shows the energy-band 
spectrum of a lattice model of Weyl fermions obtained by di¬ 
agonalizing the Hamiltonian in Eq. 56 The simple model 
of Weyl semimetal described in the previous paragraph by 
iF±(k) breaks TR, however it is also possible to realize a 
We yl syste m when TR is intact but inversion symmetry is bro- 
kei ji4 | 64 | 65 | -pjjjg implies that the system must host more than 
one flavor of pairs of Weyl fermions for the vector sum of 
kg to vanish. In the SI broken Weyl semimetal, because of 
TR symmetry there is no Berry curvature induced anomalous 
charge or thermal Hall effect in the absence of an external 
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FIG. 1. Top: Linearized band dispersion for Dirac and Weyl 
semimetals (k^ is suppressed), (a) A doubly degenerate Dirac 
semimetal (b) Transition from a Dirac semimetal to a Weyl 
semimetal (represented by a pair of Dirac cones separated by a fi¬ 
nite momentum) by breaking of time-reversal symmetry. Bottom: 
Energy band spectrum for the lattice model of Weyl fermions (kx 
suppressed) described in Eq. The two band touching points oc¬ 
cur at K± = (0,0,±7r/2). 


magnetic field. However, in the TR broken WSM, because 
of a finite Berry curvature flux through any plane interme¬ 
diate between the Weyl nodes in the momentum space, the 
ano malo us charge and thermal Hall conductivities are non- 
zerJSHEll. 

In this paper, we work with a TR broken phase of WSM 
and consider its Nemst response. Experimentally, the Nernst 
effect measures the transverse electrical response to a longi¬ 
tudinal thermal gradient in the presence of a perpendicular 
magnetic field. The Nernst effect has been used as an impor¬ 
tant experimental probe in a number of physical systems such 
as high temp eratur e cuprate superconductors^^I^, and charge 
density waveJ^El Since the TR-broken WSM has a non-zero 
anomalous Hall response (i.e., non-zero Hall efffect induced 
by the Berry curvature even in the absence of a magnetic field) 
it is expected that the anomalous Nernst response will also be 
non-zero. This is because, the anomalous Hall conductivity 
axy and the anomalous Peltier coefficient axy, which mea¬ 
sures the transverse electrical currrent in response to a lon¬ 
gitudinal temperature gradient, are related by the celebrated 

Mott relation, aab = ^ non-zero Uxy 

implies an anomalous (zero field) Nernst coefficient given by 
OLxylcTxx- In recent worlP^, however, based on a linearized 
model of a TR-broken Weyl semimetal, the anomalous Nernst 
response has been argued to be zero, because a linearized 
Weyl Hamiltonian with unbounded (or very high) ultraviolet 
cut-off of the Dirac spectrum produces daxy/dp, = 0. Here 
we show, from a lattice model of a WSM (with the lattice reg¬ 
ularization providing a physical ultra-violet cut-off to the low 
energy Dirac spectrum) that the anomalous Peltier coefficient, 
and in turn the anomalous Nernst coefficient is finite and mea¬ 


surable in a physical time reversal breaking Weyl semimetal 
such as Bii-ajSbj^. In the main part of the paper, we use the 
semi-classical Boltzmann equations in the relaxation time ap¬ 
proximation in the presence of a non-zero magnetic field and 
a Berry curvature, and derive the thermoelectric and charge 
conductivity tensors (both longitudinal and Hall) which we 
use to calculate both the conventional (i.e., magnetic field de¬ 
pendent) and topological (i.e. zero field) Nernst coefficients. 

Additionally, we also investigate the thermal conductivity 
of a WSM based on the Boltzmann equation approach. Un¬ 
like e arlier works which were based on a linearized WSM 
model,E2Ell vve employ a lattice Bloch Hamiltonian. With 
the Nernst experimental setup i.e. the temperature gradient 
VT applied perpendicular to the magnetic field B, we find 
that the transverse magneto-thermal conductivity obeys the 
Wiedemann-Franz la\v^ (i.e., the ratio of the thermal and 
electrical conductivity is the Lorenz number Lq, both with 
and without the external magnetic field). In the parallel setup 
(VT II B), however, there is an additional dependence of 
the Lorenz number, thus violating the standard Wiedemann- 
Franz law for quasiparticles in a Landau Fermi liquid, arising 
from the chiral anomaly. Our results confirm that the vio¬ 
lation of the Wiedemann-Franz law between the longitudinal 
magneto thermal and electrical conductivities^ persist in the 
physically more transparent lattice model and is not an artifact 
of the linearized low energy model. 

This paper is organized as follows: in Section II, we dis¬ 
cuss the Boltzmann semi-classical approach to calculate the 
Nernst response in a Weyl semimetal. We derive expressions 
for both longitudinal and transverse charge {aab) and Peltier 
{aab) conductivity tensors, taking into account perturbative 
electric and magnetic fields, and a finite temperature gradi¬ 
ent, for a Hamiltonian with a non-vanishing Berry curvature. 
Though the approach is general and can be applied for vari¬ 
ous configurations, we will compute our expressions relevant 
for the Nernst experimental setup. Section III concerns with 
the Nemst response in a linearly dispersing model of Dirac 
and Weyl fermions. We compute the magnetic field dependent 
Nernst response for a single Dirac node, and also for a pair of 
Weyl nodes which have a non-vanishing flux of the Berry cur¬ 
vature. Additionally Weyl fermions also exhibit an anomalous 
Nernst response even at zero magnetic field when one imposes 
a physical ultraviolet cut-off on the energy spectrum, and we 
show that this imposition gives a non-zero Peltier coefficient 
axy Section IV concerns with the lattice WSM model and its 
Nernst response. In Section V, the magneto-thermal conduc¬ 
tivity is analyzed, and the Wiedemann-Franz law is studied 
for orthogonal (VT _L B) and parallel (VT || B) setups. We 
conclude in Section VI. 


II. BOLTZMANN FORMALISM FOR NERNST RESPONSE 
IN A LATTICE WEYL SEMIMETAL 

Nernst effect measures the transverse electrical response to 
a longitudinal thermal gradient in the presence of a finite mag¬ 
netic field and absence of a charge current i.e. Ey = —-d 
dTjdx, where d is defined to be the Nernst coefficient and 











3 


—dTjdx is the temperature gradient applied along the x axis. 
The use of three conductivity tensors, a, a, and I suffices to 
relate the charge current J and thermal current Q to an ap¬ 
plied electric field and temperature gradient. We reserve the 
symbol k for thermal conductivity tensor which will be the fo¬ 
cus of Section V. We can write the following linear response 
equation 



The tensors a and a are related to each other by Onsager’s 
relation; a = Ta. In the absence of charge current (J = 0), 
we have E = d-~^a\/T. The Nemst coefficient -d can be 
derived to be 

q ddy ^xy^xx O^xx^xy 

= i-dT/dx) = ^ "aL + 

which is a function of thermoelectric tensor aab and charge 
conductivity tensor aab- We will evaluate aab and aab using 
semi-classical Boltzmann treatment in the relaxation time ap¬ 
proximation, accounting for an external magnetic field and a 
finite Berry curvature. 

A non-zero Berry curvature in a Bloch Hamiltonian acts 
like a fictitious magnetic field in the momentum spacd^, 
which substantially modifies transport properties of the sys¬ 
tem, giving rise to anomalous behavior. Anomalous transport 
due to the Berry curvature has been crucial in understand¬ 
ing intrinsic Hall and Nernst conductivity in ferromagnetic 
materialJ^mill The Berry curvature for a Bloch Hamiltonian 
H(k) is defined to be: = Cabcdk^ Ac(k), where A(k) is 

the Berry connection given by A(k) = (uk|tVk|Mk), for a 
Bloch eigenstate juk)- In the presence of Berry curvature fJk, 
the semi-cIassi cal eq uation of motion for an electron takes the 
following fortrPlEil 


and H(k) and instead write just D. In Eq.|^and Eq.|^ we have 
also defined Vk = hr^dt]^/dk to be the band-velocity. The 
second term in Eq.j^gives rise to anomalous transport perpen¬ 
dicular to the applied electric field, while the third term gives 
rise to chiral magnetic effect. The third term in Eq. [^(propor¬ 
tional to E • B) is the source of chiral anomaly, triggering neg¬ 
ative magnetoresistance. It has been shown recently that nega¬ 
tive magnetoresistance can be derived using the semi-classical 
equations of motion employing Boltzmann transporP^ Other 
recent works have also developed a modified Boltzmann equa¬ 
tion, t aking in to account Berry curvature and chiral anomaly 
effect^^SHMZMini in these works a linearized model of the 
WSM has been examined, i.e. a pair of Dirac nodes topo¬ 
logically protected by chirality quantum numbers. 

In this paper we solve the Boltzmann equation in the pres¬ 
ence of the Berry curvature and chiral anomaly terms for a 
lattice model of a WSM. The steady state Boltzmann equa¬ 
tion in the relaxation time approximation is given by 

(f • Vr-f k- Vk)/k = -^^!!^^^, (6) 

T 

where r is the scattering time, f^q is the equilibrium Fermi- 
Dirac distribution function, and /k is the distribution function 
of the system in the presence of perturbations. The scattering 
time r can in general be a function of the crystal momentum 
i.e. T = T(k), but we shall treat it as independent of momen¬ 
tum for simplicity. 

We will first consider the case when B = 0 and derive the 
longitudinal and anomalous Hall conductivities. The linear 
response relations between the charge current and the applied 
fields dictate; 

Ja = (JabEb + aab(-VhT) (7) 


f = 


1 9e(k) 
h 5k 



(3) 


where k is the crystal momentum, e(k) is the energy disper¬ 
sion, and p = hk. The first term in Eq.j^is the familiar rela¬ 
tion between semi-classical velocity f and the band energy 
dispersion e(k). The second term is the anomalous trans¬ 
verse velocity term originating from r2(k). In the presence 
of electric and magnetic fields we have the standard relation: 
p = eE -f ef X B. These two coup led equations for f and p 
can be solved together to obtairpTEII 


f = i7(B,Hk) 
p = 77(B,Hk) 


Vk + —(E X Hk) -f ■ f^k)B 

n n J 

2 

eE+J(vkxB) + ^(E-B)rJk , 

h n 


(4) 

(5) 


where i7(B,Hk) = (1 -f e(B • i7(B,Hk) is also 

the prefactor which modifies the invariant phase space volume 
dpdx —>■ D(B, Hkjdpdx, giving rise to a non-commutative 
mechanical modepTl because the Poisson brackets of coordi¬ 
nates is non-zero. For brevity of notation, we will sometimes 
omit showing the explicit dependence of Z)(B,H(k)) on B 


The charge current in the presence of an electric field and a 
temperature gradient is given bjEH 

3 = -e J [dk] (|vk +^E,x TJk) /k 

+ ^ J [c^kJJTkSk (8) 

In the above expression, [dk] = The quantity Sk = 

-/eg log /eg - ((1 - /eg) log(l - /eg)) is entropy density for 
the electron gas. The first term in Eq. [^is the current in re¬ 
sponse to an applied electric field E, also accounting for the 
transverse anomalous velocity acquired by an electron wave- 
packet due to r2(k). The second term is the anomalous re¬ 
sponse to the temperature gradient VT, which can be obtained 
using the semiclassical wavepacket methods taking into ac¬ 
count the orbital magnetization of the carriers arising from 
the finite spread of the wavefunctiorP? It can also be derived 
by first calculating the transverse heat current in response to 
an electric field and then using Onsager’s relatiorP^. The heat 
current Q takes the following form after accounting for both 
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normal and anomalous contributionJ^Ullliil 

Q = y [dk](ek - M)vk/k + J ^ ^k) Sk 

+ X J [dk]flk (^y/e, + - M)Ve,) 

- X J [dk]f2k (ln(l + + 2 Li 2 (l - A,)), 

(9) 

where Li 2 (z) is the polylogarithmic function of order 2, which 

°° k 

is generally defined as Lis(z) = X) arbitrary com- 

k=l 

plex order s, for a complex argument 1^1 < 1. The first 
term in Eq. is the standard contribution to the heat cur¬ 
rent in the absence of BetTy curvature. The second term is 
the Beny curvature mediated transverse response to electric 
field E which can be understood by the following simple ar¬ 
gument; in the presence of the Berry curvature and the elec¬ 
tric field, the electron velocity acquires the additional anoma¬ 
lous term eE x Slk- Multiplying this velocity by the entropy 
density of the electron gas, we obtain this contribution to the 
transverse heat curi'enM From Eq.[^we can write the trans¬ 
verse response of J and Q on the applied temperature gra¬ 
dient and electric field respectively as: = a^yVyT, and 

Qx = axyEy. Comparing the coefficients axy and Uxy from 
Eq.|^and Eq.[9j it is easy to note that they obey Onsager’s re¬ 
lation: Oixy = Taxy, as expected. From Eq.[^ the anomalous 
response Q on an applied temperature gradient can be written 
as: Qx = Ixy^yT. The quantity Ixy in Eq.[T]can be calculated 
as: Ixy = —kgTc 2 /fi, wherdST 


J, and comparing the resulting expression with Eq. the lon¬ 
gitudinal components of the conductivity tensors dab and aab 
can be easily read to be 


o-xx = e^y [dkJu^T (13) 

axx = -^ J [dkjvlrie-fi) (14) 

where Vx = h~^d€\i/dkx- The transverse components are; 

CTyx = y y [dkjflfeq, (15) 

ayx = ^ I [dkJJTkSk, (16) 


which are purely anomalous because we have assumed E and 
VT are applied along the x direction, and there is no magnetic 
field. 

We now discuss the case of a finite magnetic field. We con¬ 
sider a particular configuration relevant for the experiments 
measuring Nemst coefficient i.e. VT = WxTx, B = Bz, 
and E = 0, although the approach will work for other con¬ 
figurations also, like the parallel setup discussed in Section V. 
The Boltzmann equation (Eq.j^takes the following form after 
making substitutions for f and p from Eq.j^an d|g. 


VxTVxT 


e-M 

T 



/k - feq 

D{B,nk)T 


eB f d d \ 


/k 


(17) 


oo 

c„= j [dk]n, J de{(3er^ 

€-/i 


( 10 ) 


The following Ansatz is chosen for the distribution /k which 
also accounts for correction factor (A) due to a finite magnetic 
field. 


The e nergy integral in Eq. fTOlreduces to the following for n = 

jma ' 


/k = fey -[ DrVxVxT^-dJ^ _ V . A U - 


T 


de 


(18) 


de(/3e) 


de 


— ~^feq + / 3 ^(e — df feq 


e-fi 


(ln(l + + 2Li2(l - feq)) (11) 


Eq. 11 and Eq. 10 combined with Ixy = —kgTc 2 /h give the 


last two VT dependent terms in Eq.[^ 

Keeping only linear order dependence on the applied field 
E and VT, the following Ansatz is assumed for the distri¬ 
bution function, which is a solution to the steady state 
Boltzmann equation (Eq|^ : 


fu= feq +t{- 


dfeq 

de 


Vk 


—eE 


e-M 

T 


(-VT) V(2) 


We will further assume that the electric field and temperature 
gradient have non-zero components only along the x direc¬ 
tion. Substituting for /k from Eq. 12 in Eq.|^for the current 


The Boltzmann equation (Eq. 17 1 thus becomes 
d 


eB 


d 




V A 


(19) 


Imposing the condition that this equation must be valid for all 
values of v, we find that A^ = 0, and the equation can be 
simplified to; 


cB^xT^-tJ^Dt 


T 

- '^xEx I 


Vx 


V, 


V 


•’xy 


m 

eB 1 

-1- 

Dt 


m 


■xy 


-VyAy 


I _ '^xAy 

( 20 ) 


'^XX 

eB 1 


m 


xy 


Dr 


In order to solve the above equation, we introduce complex 
variables V = ivy, and A = — iAy, and rewrite the 
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equation in the following manner 
Re 


eBrDVrT^-^V 


= Re 


T 

VK\—-^\ + 

\ Dt / 


+ 


eBVA* 


m 


xy 


( 21 ) 


where = h '^d'^E{k)/dkidkj is the inverse band-mass 
tensor and Re( 2 ;) stands for the real part of z. Eq. |^can be 
solved for A: 


Aa; = eBTD{B,^l\^)V xT- 



t'x '’y 


eBv^ 

1 eBvx 

Vx 



Vx _|_ 

Vy 


eBvx 

- k 

rrixy m^x 



rrixy 

Dt 


rrixx 

TTlxy 


rrixx 

T 



eBvy 

1 eBv„, 

Vx 

2 

eBvx 

eBvy 

Vy 




rrixx 

rrixy 

Dt 



ruxx 

rrixy 


Dt 


eBv^j 


Dt 


( 22 ) 


Ay = eBTD{B, r2k)V2;T- 


T 





eBvy eBvy: Vy 


l^X 


eBvy eBvx -Ux 

k 

TTLxy rrixx 


rrixy ' rrixx Dt 


rrixx rrixy 


mx:x rn,cy Dt 


eBvj. 


Dt 


Dt 


(23) 


For convenience of notation, we rewrite A^; and Aj, as; Aj = 
tV xT^^Ci, incorporating into Ci the remaining factors apart 


from tWxT^j^ of Eq. 


22 


23 


Using Eq. 


18 


and the 


and Eq. 

results for A^: and Ay, w^can now explicitly writelHe distri¬ 
bution function /j. as: 


/k = feq-{ tV,T 


A- k- (df, 


T 


I eq 


de 


((Ca; B)Vx CyVy) 

(24) 


The expression for the charge current J, in the pr esenc e of B 
and flk, is also modified by the factor D(B, as we 

pointed out earlier that D(B, Uk) = (1 + e(B • Uk)/^)”^ 


IS 


the multiplicative factor which alters the phase space volume 
locally. 


J = —e / [dk]Z? ^f/ 


kseVT 


[dk]f2kSk (25) 


Substituting Eq. 24 in Eq. and again comparing with the 


atj can be solved to: 


linear response relations in Eq. the thermoelectric tensor 

dfeq 


0!.xx — C 


[dk]vx 


T 


de 


(cx - D) 


ayx=e J [dk]{vlcy -f (cx - D)vxVy) 

[dkJO^Sk 


kse 

h . 


(26) 
dfeq 

de 
(27) 


The temperature dependence of the zero-field anomalous con¬ 


tribution in Eq. 27 is hidden in the entropy density Sk of the 
electron gas. Similarly, the electrical conductivity compo- 


I 

nents (transverse and longitudinal) are obtained to be: 

dfeq' 


[dk]r 


de 


(cx - D) 


(28) 


J [dmvlcy + VxVyiCx - D))t 

f [dk]n^feq 


(29) 


As a good check of our calculation we also recover the results 
for aab and a^b found in Ref. |45] where the tensorial nature 
of rriij is ignored and m = h^/v^. The transverse compo¬ 


nents, i.e. Eq. 27 and are a sum of two terms: the first 
term captures the effect of a finite B which is further modi¬ 
fied by the Berry curvature flk, due to the factors Cy, Cx and 
D which are non-trivial functions of the Berry curvature. The 
first terms in Eq.j^and Eq. [^depend on the scattering time 
T, and we call these as ‘modified’ i?-dependent Hall conduc¬ 
tivities (because they are modified due to the Berry curvature). 
The Berry curvature also alters the expressions for longitu¬ 


dinal conductivities given in Eq. 26 and Eq. 28 because of 
the factor Cx — D. In the limits when flk —0, the factor 
Cy WT for a quadratic band dispersion, upto linear order 
in B, (where w = eB/m is the cyclotron frequency). In the 
same limit, the factor Cx — D ^ —1 upto zeroth order in B, 
thus yielding the standard expression for (Jxx and Uxx given 
in Eq. 13 and 14 In contrast, the second term in Eq. 27 


and 29 which is Berry curvature dependent persists in the 


absence of a magnetic field, and is a purely anomalous contri¬ 
bution. We shall roughly examine the limit in which the factor 
i9(B,Uk) = (1-f e(B • Uk)/S)“^ significantly deviates from 
1. Defining k = 27rK/a, where a is the lattice constant, and 
K is dimensionless, e(B • = a^iln/lg, where Is is 

the magnetic length and Uk is dimensionless. For a magnetic 
field of IT, a - 21, e(B • n^)/h ~ IQ-^VIk. The Berry 
curvature for a single linearly dispersing Weyl node centered 
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at the origin is JIk = K/4|Kp, thus D{B, 17k) ~ 1 for low 
magnetic helds and away from band touching point (Kq = 0) 
in the momentum space (i.e. approximately |K| ^ 0.01 in 
this case). This is expected, as qualitatively one understands 
that the effect of Berry curvature peaks when the energy band 
gap Eg —> 0. When the effects of Berry curvature can be ne¬ 
glected, the following standard expressions are derived from 
Eq.[^and|^for Hall conductivities, keeping terms only upto 
linear order in B: 


^ xy 




VxVyd'^e 

dkxdky 

(30) 


v'^d'^e VxVyd‘^e\ 

dkl dkxdky J 

(31) 

In Section III, we will use the formula obtained for aab and 
Uab to calculate the Nernst coefficient in Eq. hrst analyti¬ 
cally for a simple Dirac and Weyl linearized Hamiltonian, and 
then numerically in Section IV, for a lattice model of Weyl 
fermions. 




Oi-xv — 


Th 


J [rfk](e-p) 


III. NERNST RESPONSE IN LINEARIZED MODEL 
DIRAC AND WEYL SYSTEMS 

In this section, we will concern ourselves with the Nernst 
response of a linearized spectrum of Dirac and Weyl sys¬ 
tems. We examine the magnetic field dependent transverse 
conductivities (axy and axy) for a linearly dispersing Dirac 
node, which are analytically tractable using the Boltzmann 
approach. We repeat the procedure for a pair of Weyl nodes 
taking into consideration the Berry flux modification of the 
normal B-dependent conductivities. 


A. Nernst effect in a linearized Dirac Hamiltonian 


As a warm up, we discuss the Nernst response of a single 
Dirac cone, with linear dispersion Ck = where k = 

^kx -f -I- k^. The density of states for a single Dirac/Weyl 
node with unbounded linear dispersion (taking into account 
spin degeneracy) is given by 


p[E) 


1 E^ 

TT^ 


(32) 


The density of states vanishes at the Dirac node, which gives 
rise to many unusual properties. Eq. for the longitudi¬ 
nal conductivity, without the Berry curvature term, reduces 
to Eq. [T^ which can be employed to analytically deduce the 
zero temperature conductivity for a Dirac Hamiltonian to be 


9 2 


( 33 ) 


where we have assumed the scattering time r to be a phe¬ 
nomenological parameter independent of energy or momen¬ 
tum. Eor transverse magneto-conductivity, for a weak mag¬ 
netic held B = Bz, we derive the Hall conductivity axy for a 
linearized Dirac Hamiltonian using Eq.|^to be 


e^Br'^ vfP 


(34) 


We note that {axyjaxx) = wr = eBrlm, where m is the 
band mass near the Eermi surface given by m = Hp/vp. At 
low temperatures, the thermoelectric tensor aab is related to 
the derivative of aab via the Mott relatiorl^ 



k%T dagb 
e dp, 


(35) 


Combining the Mott relation with Eq.|^ the Nenrst coefficient 
becomes (when axx ^ cTxy) 




dp 


3Hp'^ 


(36) 


36 


IS 


where Qh — ^^xyj^xx is the Hall angle. However Eq. 
valid only when p ^ 0, and d does not diverge for p = 0 as 
we shall see shortly. 

The scattering time in Boltzmann conductivity is sensitive 
to the type of impurities in the system. Eor neutral short range 
or point-like impurities, the scattering time Tg is given byP^ 


J_ ^ UsV^kl 
Ts 2lFh?VF ’ 

where Ug is the density of the impurities, Vq is the strength of 
the impurity potential, and kp is, the Eermi wave-vector. Con¬ 
sidering Thomas-Eermi (TE) screening, the scattering time for 
long-range ionic impurities at zero temperature is given bjEH 

— = Aira'^na^Itiqo) (38) 


In the above expression, ric is the density of charged im¬ 
purities, a = /kHvf is the hne-structure constant, qq = 
qpF/'^kF, where qpF is the Thomas-Eermi wave-vector, and 
It{x) = -f l/2)log(l + 1/a;^) — 1. The total scattering 
time is given by Matthiesen’s rule 


1 1 1 

- = — -f — 

T Ts Tc 


(39) 


Eor a linear Dirac Hamiltonian, kp = p/hvp, therefore the 
scattering time expression take the following form: 

_ nc4TTa'^v^hIt{qo) 

Tc p^ 

_ usVqp'^ 

Ts SttH^Vf 


(40) 

(41) 


Erom Matthiesen’s rule, it is evident that the shorter time scat¬ 
tering process (t^ or Tc) will dominate the carrier transport, 
therefore near p = 0, the scattering from ionic impurities will 
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primarily determine the conductivity, and for /r ^ 0, it is 
scattering from the neutral point-like impurities that govern 
charge transport. For an arbitrary value of /i, the following 
expression for scattering time r can be written, using expres¬ 


sions in Eq. 39 40 and 41 


1 1 / 1 . d 

- = — ^(1 + a;/ 

r To Vfi 


(42) 


where tq and x are constants depending on the coefficients of 
and l//r^ in Eq. 40 and 41 whose exact form is lengthy 


and not illuminating for our discussion. Using this expression 
for the total scattering time t in Eq. we obtain: 




Ott^ VFh^(l + xfi‘^) 


Similarly, from Eq. 34 we have 

eBvpTofJ,^ 


^ XV — 


Ott^ h'^(i -\- 


The Hall angle near /i = 0 reduces to: 

0/r = eBfivpTo/H, 


(43) 


(44) 


(45) 


and thus from Eq. 36 the Nemst coefficient at /i = 0 is 
given by: 


t?n = - 


TT^ k'^T eBvpTQ 

3 e h 


(46) 


Alternatively, the same conclusion can be reached by using the 
equation for the Nernst coefficient i.e. Eq.|^ and Eq.|^ 43 
and|^ by substituting the exact expressions instead of using 
the Mott relation. Ear away from the Dirac point at /r = 0, the 
charge conductivities are 


^XX 


^ xy 


'To 





eBvpTQ 


(47) 

(48) 


The Hall angle in this case no longer varies linearly with the 
Eermi energy (as in the case near /r = 0), but is instead given 
by Qp = eBv'^To/hx^. The Nernst coefficient for /r ^ 0 
becomes 


t? 


ir^kpT eBvpTo 


(49) 


which approaches zero as fi is increased asymptotically, as ex¬ 
pected from Eermi liquid theory, where the Nemst coefficient 
vanishes because of Sondheimers cancellation.!^ The sign of 
the Nernst coefficient does not change with the sign of p. and 
is thus an even function of p. This is because both a^yOCxx and 
o'xxO/.xy, which appear in the numerator of the expression for 
-d in Eq.j^do not depend on sgn(p). The plot in Eigure [^dis¬ 
plays Nernst coefficient for a linear Dirac node obtained using 
Eq.[ 


and 35 We have provided a physical ultra-violet 
cut-off to the low energy spectrum at e = ±3.5f with t = 0.1 




FIG. 2. Left: Plot of the Nemst coefficient as a function of p/f 
for a single linearized Dirac node whose spectrum is bounded at ± 
3.5f. The energy parameter t is chosen to be f = 0.1 eV. Here i9o 
indicates the value of the Nemst coefficient at p = 0. Right: Chosen 
scattering time for the calculation in ps as given by Eq.|42| 


eU, vp 10® m/s T = 40K, and B — IT. The chosen 
scattering time t used for this calculation, and its scaling with 
p is also shown in Eigure [^ A regularized lattice model will 
however smoothly bound the dispersion in the Brillouin zone. 
Section IV will be devoted to evaluating the Nernst response 
of a WSM Hamiltonian defined on a lattice. 


B. Nernst effect for a pair of linearized Weyl fermions 


A single Dirac node can be visualized as two Weyl nodes, 
which are topologically protected by chirality quantum num¬ 
ber of opposite sign, coinciding with each other in energy- 
momentum space. Therefore the net flux of the Berry cur¬ 
vature, and henceforth the net chirality vanishes for a single 
Dirac node, resulting in a zero anomalous response for both 
the charge and thermoelectric conductivity. As a result, no 
anomalous Nernst response is also expected in a linear Dirac 
Hamiltonian. Using an external perturbation, a single band¬ 
touching point in a Dirac cone can be shifted into a pair of 
isolated Weyl points possessing opposite chirality quantum 
numbers. The external perturbation must break either time- 
reversal symmetry or inversion symmetry, and also lifts-up the 
degeneracy of the Dirac spinor. Our discussion will be cen¬ 
tered upon the assumption that time-reversal symmetry is vio¬ 
lated, which can be achieved using magnetic field as a pertur¬ 
bation. Also the Weyl points are assumed to occur at the same 
energy, thus there is no chiral chemical potential. One can 
also construct inversion asymmetric and TR invariant models 
of a WSM, but they are not of interest to us here because the 
anomalous Hall and Nernst response vanishes as the vector 
sum ko of the node separation becomes zero. Eigure]^ shows 
Berry curvature plot of a Weyl semimetal in the kz = 0 plane, 
where we assumed the node separation is kg = (0.5, 0.5,0). 
At the origin is a Weyl node with chirality quantum number 
H-1 which acts as a source of Berry flux. At ko we have another 
node with chirality quantum number -1 which acts a sink of 
Berry flux. 

Let us first concern ourselves with the normal contribution 
to the Nernst effect i.e. due to an external magnetic field. As 
we pointed out earlier, this contribution is further modified 
due to effects of the Berry curvature, which are encoded in the 
factors Cx, Cy and i9(B, Uk), expressed in Eq. 26 ■ 29 Adding 
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FIG. 3. (color online) Plot of Berry curvature of aWeyl semimetal in 
the ^2 = 0 plane, where the node separation is ko = (0.5, 0.5, 0). At 
the origin is a Weyl node of positive chirality, which acts as a source 
of Berry flux indicated by outgoing arrows, which flow towards ko, 
which acts as a sink of Berry flux indicated by incoming arrows. 


up the contribution from both the Weyl nodes and keeping 
terms only upto linear order in the magnetic field B = Bz 
(which is justihed in the limit lot <C 1), the zero tempera¬ 
ture longitudinal electrical conductivity for the Weyl system 
becomes 

I - M) 

r sin^ ed9 
Jo ~ h'^e'^VpB'^ cos^ 9 


For B ~ IT, vp ^ lO^m/s, (which is the typical value for 
vp in a WSIViS]), the factor fi^e^VpB^ ~ 10 Thus away 
from fi = 0, the expression reduces to Eq. 33 for a Dirac 
node, except upto an overall multiplicative factor of 2 (for two 
nodes). At p = 0, where the Fermi surface just reduces to a 
pair of Weyl points, a^x = 0- From Eq. 29 we can calculate 
the charge Hall conductivity axy for a WSM (again upto linear 
order in B): 


^yx 


= 2e2 



VyUiT'^{l + (eBilz/h)^) 

(1 - {eB^.jhYY 


|J(ek - 9) 


sin^ 9fj,^{fi^ + fi^e^VpB^ cos^ 9) 
Ap'^h?vp Jq — h?e^VpB^ cos^ 9)^ 

(51) 


Further away from p, = 0, the above expression also reduces 
to Eq. 1^ for Dirac node, after making the substitution for 
the cyclotron frequency oj = eVpB/hii. At p = 0, again 
c^xy = 0- We can therefore conclude that away from the band 
touching point (which line-up with the hne tuning of chemi¬ 
cal potential p = 0), and for typical values of the Fermi ve¬ 
locity and weak magnetic field (such that the semiclassical 
Boltzmann approach is still valid) the deviation of the nor¬ 
mal Nernst response due to the Berry curvature is negligible. 
In this limit the normal contribution to the Nernst response 
roughly reduces to the sum of individual contributions from 
the two Dirac nodes. We have also verified this conclusion 


through explicit numerical integration, even at hnite tempera¬ 
tures, using Eq. [26l - [29| 

A Weyl system also exhibits anomalous Nernst response 
even at zero-field and thus the total Nernst signal must arise 
from both contributions. The anomalous Hall conductivity 
a^y for a time-reversal broken Weyl semimetal is non-zero 
and varies linearly with the node separation fcoP^which can be 


obtained by integrating Eq. 29 with the corTect regularization 


that is consistent with the broken symmetries in the presence 
of fco.® 


h 27r2 


(52) 


The result in Eq. 


52 


suggests that the remains unaltered 


with temperature or for a hnite p. This result is strictly valid 
only for an unbounded linear dispersion of Dirac fermions. 


Mott formula (Eq. 35|l then suggests that a^y = 0. If an 


upper physical cutoff on the energy of a Dirac node is im¬ 
posed then a^y is non-zero, because the contributions from 
the partially hlled states generically will remain hnite. Fig¬ 
ure 0 shows the plot for a^y as function of node separation 
ko for a linear WSM with an upper energy cut-off, and also 
the plot for a^y as a function of upper energy cutoff, obtained 


using Eq. 27 and 29 Lowering the cutoff results in a hnite 


non-zero axy, thus the anomalous Nernst response is also ex¬ 
pected to be non-zero. 

It is not entirely evident from Eq. that the anomalous 
Nernst response will vanish for a linearized Weyl Hamilto¬ 
nian with an unbounded dispersion when a^y = 0, because of 
the non-zero factor axxC^xy the numerator of Eq. How¬ 
ever, using the results from the previous subsection we have 


from Eq. 


43 


O5 ^xx — rr. 


^ "xx9^ 


and from Mott 


that near p 

relation: Uxx = ctxx9^’ thus in the vicinity of the Dirac point, 
the anomalous Nernst coefficient becomes: 


■9 = 






(53) 


Two points can be noted from the above equation: the anoma¬ 
lous Nernst coefficient vanishes at the Dirac point when p = 
0, and the Nernst coefficient is an odd function of p. The 
evaluation of total Nernst signal, will have contributions from 
both normal and anomalous Hall conductivities, and therefore 
an asymmetric behavior of the total Nernst coefficient about 
p = 0 is expected. 


IV. NERNST EFFECT IN A LATTICE WEYL 
HAMILTONIAN 


The novel semimetallic state with Weyl-like fermionic exci¬ 
tations has been recently realized in a series of experiments. In 
the inversion asymmetric crystalline compound TaAs, exper¬ 
iments have claimed to host topologically non-trivial Fermi 
arcs and Weyl cone JSHUl Another pathway of observing Weyl 
fermions in condensed matter is to break TR symmetry by ap¬ 
plying a magnetic field in a 3D Dirac semimetal and split each 
Dirac cone into a pair of Weyl nodes. 3D materials NaaBi, 
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we consider a prototype lattice Hamiltonian for 
Weyl fermions described by in Eq. 56 because it 

serves our purpose of discussing the Nernst effect. 


HiattCk.) = t{sin{kxa)aj; + sin{kyb)ay + cos{kzc)az) 

+ m{2 — cos{kxa) — cos{kyb))a^ = ■ cr, 

(56) 


FIG. 4. Left'. Plot of anomalous a^y (in the units of fikse/K) for a 
bounded linear WSM with two bounded linear Dirac nodes vs. the 
upper energy cutoff Ec (in the units of lOmeV), for fj, = —0.7f and 
fco = 0.25. As the upper energy cutoff for each bounded Dirac node 
increases, a^y decreases eventually becoming close to zero. Right. 
Anomalous Hall conductivity a^y as a function of the node separa¬ 
tion ko for a bounded linear WSM, where is the Hall conductiv¬ 
ity when ko = 0.15. 


Cd.-j As 2 , and Bii-j^Sba; with a; ~ 3 — 4% were recently pro- 
poseci^32iim®£^cl also realized experimentally to be Dirac 
semimetal^2lEEI] xhis has paved the way to realizing a TR 
broken WSM phase, w hich has b een so far verified by a few 
experimental signature^^HlIEIlIl 

One can construct an effective Hamiltonian Hjj (k) describ¬ 
ing electron dynamics near the symmetry points in momentum 
space of a 3D Dirac semimetal 

Hoik) =VF{'k-a)T:^+m{k)Ty„ (54) 


Hiattik) supports a pair of Weyl fermions located at 
(0, 0, ±7r/2c),thus ko = (0,0, tt/c) when m > t/2. This lat¬ 
tice model of Weyl fermions described in Eq. can mimic 
the experimentally relevant TR breaking WSM model de¬ 
scribed in Eq. and the paragraph below it, when we con¬ 
sider an external magnetic field B = (0, OjTtuf/p), as the 
node separation is fixed from Eq. This is however suffi¬ 
cient to discuss Nernst effect in the lattice WSM. The energy- 
band spectrum is shown in Eigure[2 obtained from diagonal- 
ising the Hamiltonian in Eq. The anomalous Hall con¬ 
ductivity cr^ for this Hamiltonian is given by Eq. 


15 


at finite 

temperature and chemical potential, which is plotted in Eig- 
ure^ At zero temperature and at /r = 0, a^y = ^ jh. The 
component of BetTy curvature vector Jlk, for Hiatt (k) is 
given by 


n 


k,n,i 


Nk- 

{-iTe^yi - 


dkj dki 


4|Nk|3 


(57) 


where a and r are the vectors of Pauli spin matrices acting 
on the spin and pseudo-spin degrees of freedom respectively. 
The mass term m(k) can be tuned to vanish at isolated points 
in the Brillouin zone by choosing a specific form of m(k) as a 
function of crystal momentum k. The mass term m(k) can be 
chosen to be TO(k) = TO+pcos(k), for |k| ^ 1, which can be 
realized physically. Eor instance, at a; ~ 3 — 4% in Bii_ 2 ;Sb 2 . 
the mass term can be tuned to zero at particular K values. This 
yields a degenerate linearly dispersing 4-component Dirac 
fermion described by Eq. 54 with to(K) = 0. A perturba¬ 


tive Zeeman field coupled with the spin degree of freedom 
described by Hz = gBzCz, can now lift the degeneracy of 
the Dirac fermion, with the energy spectrvim given bjHH 


E{k) = ±^vjz{kl + kl) -F {gBz ± vpkzY, (55) 

where the band-touching point is shifted from E point at 
K = (0,0,0) to ±Ko = -ii{Q,Q,gB /vf), where g is the 
Lande g-factor. The direction along which the two nodes ap¬ 
pear is along the direction of the applied magnetic field. The 
separation /cq between the two Weyl nodes is magnetic field 
dependent and is given by fco = 2gB/vF- Now the Hamilto¬ 
nian around each Kq is that of a linearized Weyl fermion. 

To discuss the Nernst response in a physical Weyl system, it 
is advantageous to consider a lattice model of Weyl fermions 
with the lattice regularization providing a physical ultra-violet 
smooth cut-off to the low energy Dirac spectrum, because 
the linearized continuum theory for calculating the anomalous 
Hall current at a finite density turns out to be insufficient. Eor 
simplicity, instead of considering the band Hamiltonian for 


In Eq. 57 n stands for the band index. The anomalous Peltier 
coefficient , can be calculated using Eq. 


16 


and is non¬ 
zero, as shown in Eigure N because a^y is no longer a constant 
function of g. Eigure p also shows a^y calculated using the 
Mott relation. This feature of non-zero was absent in the 
linearized model of a WSM. Eigure m shows the total Nernst 
response d obtained by numerical calculation using Eq.[^p^ 
|29]for this lattice model. We observe that the normal Nernst 
coefficient is non-zero around p, = 0, and is an even function 
of p, which consistent with the findings of Sec. III. In Sec. 
Ill, we also pointed out that the anomalous Nernst coefficient 
vanishes at p = 0 and is an odd function of p, and thus the 
total Nernst response is expected to show an asymmetric be¬ 
havior about p = 0. These conclusions do no change for a 
lattice model, and a slight asymmetry can be observed in the 
plot of Nernst coefficient in Eigure]^ 


V. MAGNETO-THERMAL CONDUCTIVITY FOR A WEYL 
SEMIMETAL 


When a temperature gradient VT is applied across a sam¬ 
ple, a charge current J is developed. In the absence of a charge 
current, from Eq. we must have E = a~^6iWT. Using 
this expression for E in the formula for the thermal current 
Q given in Eq. [T] we can write the following linear response 
relation 

Q = {Taa-'^a - i)VT = k{-VT) (58) 
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the quasiparticle description of electronic states remains valid, 


‘^=LoT 

O'ij 


( 61 ) 


Eq. 61 states the Wiedemann-Franz law, where Lq is the 


Lorenz number (Lq = 7r^fc^/3e^). 

In the absence of a magnetic field, the B-dependent con¬ 
tribution to axy is zero, and k^x = Ixx — ct^xxl'^xx will be 
given by the standard expression of the longitudinal thermal 
conductivity: 


Ixx = J [dk]vl 


T 


dfe. 


de 


(62) 


To discuss magneto-thermal conductivity, we first examine the 
case which is also relevant to the Nernst experimental setup 
discussed in Section II, III and IV i.e. VT = {dT/dx)x and 
B = Bz. The expressions for charge and thermoelectric con¬ 
ductivities a and a have been already obtained in Eq. 26p9 


From Eq.j^and Eq.[^ we can read the conductivity tensor I 


Ixx = J [dk]vl 




dfeq 
' de 


(cx - D) (63) 


FIG. 5. (color online) Results for the lattice WSM described by 
Eq. Top Panel (left)'. In red is the anomalous contribution to 
e^xyloi^y obtained using Eq.flm, and in blue is the anomalous con¬ 
tribution to Uxy using Mott relation (Eq.|35[l displaying a reasonable 
agreement; q™ is the value at p, — t. Top panel (right): Anoma¬ 
lous Hall conductivity as a function of p/t (again a^y is the value at 
p = 0)) obtained using Eq. 15 Middle panel (left): Normal Nernst 
response i9^/i9o as a function of p, where do is the Nernst coeffi¬ 
cient at /i = 0. Middle panel (right): Anomalous Nernst response 
d^/dm as a function of p, where dm is the anomalous Nernst co¬ 
efficient at p — 0.3t. Bottom panel (left): Total Nernst coefficient 
d'^ j dm including both normal and anomalous response as a func¬ 
tion of pjt, where dm is the Nernst coefficient at p = O.lt. A slight 
asymmetry for the total Nernst signal about p = Q can be noted. Bot¬ 
tom panel (right): Chosen scattering time in ps as a function of p, as 
given by Eq. The chosen parameters values for this calculation 
are; T = 2QK, t = 0.1 eV, m = 0.6t. 


k is the thermal conductivity tensor whose longitudinal and 
Hall components can be written explicitly as: 


^Xx{,tXxx ^xy^ ‘^tJxyCtxxCtx 


(59) 


'*^xy 


— ^xy 


'xy\ 


(o^ly - aL) + 


■' xx^xx^xy 


xy 


(60) 


Usually I is identified with k, however the second term in 


Eq. 59 and Eq.|^is still non-zero, though small compared to 
1. Wiedemann-Franz law states that the ratio of thermal con¬ 
ductivity K and electrical conductivity a for a metallic state is 
proportional to temperature. The law holds for a generic sys¬ 
tem as long as it can be termed as Landau Fermi liquid where 





- X J [dk]nk(ln(l + + 2 Li 2 (l - fey)) 

(64) 

The B-dependent longitudinal magneto-thermal conductivity 
Ixx is further modified from Eq. |^by the factor of (c^ — D) 
which is a function of the Berry curvature. The first term of 
the transverse magneto-thermal conductivity lyx in Eq. [64|is 
the standard B-dependent contribution. The second and the 
third terms give the zero magnetic-field anomalous thermal 
conductivity. The Wiedemann-Franz law given in Eq. [^re¬ 
mains valid for Kxx and Kxy as shown in Figure]^ 

A more interesting scenario occurs when VT = [dT / dz)z, 
E = 0, B = Bz i.e when the applied temperature gradient 
is parallel to the magnetic field. Using Eq. |^and Eq. the 
steady state Boltzmann equation (Eq.|^ becomes 



The following Ansatz is chosen for the distribution function 
/k, which is a solution of Eq.|^ 

/k - fey = -Dt^-^V,T {v, + |t(v • f^k)) 

+ (-f^)v.A (66, 
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fl/t 


FIG. 6. (color online) Top panel (left)'. Longitudinal magneto- 
thermal conductivity K:cx for the WSM lattice model as a function 
of p./t, where is the value at p = t. Top panel (right): Normal 
B-dependent contribution to where is the value at p = t. 
Bottom panel (left): Anomalous zero magnetic-field contribution to 
where K%y is the value at p — 0. In all the three figures, we 
have plotted in red the thermal conductivity obtained using Eq. 
and Eq.|^ and in blue using Wiedemann-Franz law Eq.|^(all three 
plots are for transverse setup i.e. VT _L B). Bottom panel (right): 
Plot of AL{B)/Lo = (L(B) /Lq — 1) as a function of applied mag¬ 
netic field B (here Bq = VFko/2g) when VT || B showing an addi¬ 
tional B^ dependence of the Lorenz number arising from the chiral 
anomaly term (f2 ■ v). 


tudinal conductivities are obtained to belSEH 

eB 


T 


D[dk] 


O 


Izz = I D[dk] ( Vz + • V 


(e - fi) - 




T 


dfeq 

de 
(70) 


dfe. 


de 
(71) 


In a similar fashion one can also calculate the charge conduc¬ 
tivity azz- 


ctzz = e" / D[dk] ( Vz + r 


eq 


(72) 


The results in Eq. [70]-Eq. [7^ are of interest because of the a 
dependence arising from the chiral-anomaly term 17 • v. 
Note that this term does not arise in the transverse setup i.e. 
when VT and B orthogonal to each other, and is thus linked 
to the topological E • B term arising in axion-electrodynamics 
of WSM. We can write the following simple relations for crzz 
and Kzr- 


cTzz = 0-0 + aB^ 

Kzz = K 0 + PB'^, 


(73) 

(74) 


The correction factor A in the Ansatz for /k is introduced to 
account for a perturbative magnetic field B. Substituting for 
/k given in Eq.[6^into the Boltzmann equation (Eq.[65|, and 
imposing the condition that the equation should remain valid 
for all values of v, we find that = 0. Introducing V = 
Vx + iVy and A = — iAy, the Boltzmann equation can be 

rewritten in the following form 

Re I 

eBDri^e — p)AJzT f h 
rriyz 

+ 

hlTLyx Dt 


-LRe ( V 
= -Re 
where. 


Th? 

-ieBVA eBV*A 


17^ 
ttlxx 
17 T, 


17,, 


•’xy 


S = 


m 


xy 


17,, 


ruy 


17z 

ttlxz 

172 

rriyz 


(67) 

( 68 ) 

(69) 


The factors A^; and Ay can be straightforwardly evaluated 
from the real and imaginary parts of the complex vector 
A = Aj, + iAy, which is a solution of Eq. 67 Now sub¬ 


stituting for the distribution function /k in Eq. 25 it is then 


where erg and ko are the longitudinal conductivities for B = 0. 
The coefficients a and (3 (not to be confused with (3 = 
l/ksT) account for the B^ dependence of azz and Kzz re¬ 
spectively. They depend on the band-structure of the Hamil- 
tonain and can be obtained by the k— space integrals defined 
in Eq. ■ 71 The Lorenz number L in the Wiedemann-Eranz 
law given in Eq. 61 will be B-dependent, and can be written 


L{B) = 


AL{B) 


= Lo + AL{B) = 


kq + (3B^ 


TcJzz 
(/3tTo - QfKo)B^ 
Tal 


T(ao + aB2) 

(/3cro - aLoTao)B‘^ 


Tal 


(75) 


(76) 


AL{B) gives the B^ enhancement of the Lorenz number from 
its standard value Lq. Eigure|^displays the quadratic behavior 
of the Lorenz number L{B) with the magnetic field in the par¬ 
allel setup for the lattice WSM Hamiltonian given in Eq.[56l 
thus showing a violation of the Wiedemann-Franz law. It is 
also worthwhile to point out the sign of AT in Eq. 76 which 


will depend on the details of the band structure of the Hamil¬ 
tonian. The Lorenz number, which is the ratio of thermal to 
electrical conductivity will increase (decrease) from its stan¬ 
dard value if the B^ coefficient in the expression for thermal 
conductivity is greater (lesser) than electrical conductivity In 
the present case, the sign of AT was found to be positive. 
Similar conclusions on the sign of AT were obtained in pre¬ 


possible to deduce the conductivity tensor a and 1. The longi- 


vious won 


m 
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VI. CONCLUSIONS 

In this work we have studied the Nernst response of a time- 
reversal broken Weyl semimetal. As a consequence of non¬ 
zero anomalous Hall response in a Weyl system with bro¬ 
ken time-reversal symmetry, it is generally expected that an 
anomalous Nernst conductivity is also observed. This is be¬ 
cause generally the Peltier coefficient which is related to the 
first derivative of the charge conductivity with respect to the 
chemical potential should not vanish. However, a linearized 
Weyl fermionic system was found to have its anomalous Hall 
conductivity independent of chemical potential and tempera¬ 
ture. Previous studie^have therefore argued that the anoma¬ 
lous Peltier coefficient and the anomalous Nernst response for 
a system of Weyl fermions should be zero. We show this by 
considering a physical description of a WSM which is cut¬ 
off at higher energies by either considering a bounded lin¬ 
earized Weyl Hamiltonian, or a lattice regularization provid¬ 
ing a smooth physical ultra-violet cut-off in a lattice model of 
Weyl fermions. This produces a non zero Peltier coefficient 
and thus a non-vanishing Nernst response measurable experi¬ 
mentally. 

Starting with the semi-classical Boltzmann approach to lin¬ 
ear transport in a system, we hrst derived the expressions for 
charge and thermal conductivities in the presence of a pertur¬ 
bative magnetic field and a temperature gradient orthogonal 
to each other, for a generic band Hamiltonian which has a 
non-trivial Berry curvature. The longitudinal conductivity is 
modified from its standard expression because of Berry curva¬ 
ture effects. The B-dependent transverse conductivity also is 


modified by Berry curvature. Additionally, the transverse con¬ 
ductivity also comprises of a purely anomalous contribution 
even at zero B-field due to the Berry curvature. Thus the to¬ 
tal contribution to the Nernst signal comprises of two parts: a 
B-dependent response, and a purely anomalous response. We 
derived analytic expressions for the Nernst coefficient in a lin¬ 
earized Dirac and Weyl Hamiltonian, and have also computed 
the total Nernst response for a lattice model of Weyl fermions 
numerically. We also pointed out that B-dependent normal 
Nernst signal is an even function of the chemical potential, 
but the anomalous Nernst coefficient is an odd function. As a 
result one would expect a slight asymmetry in the total Nernst 
response p = 0, which is evident from our numerical studies. 

Additionally, we also examined the magneto-thermal con¬ 
ductivity of a WSM, and hnd that for orthogonal experimen¬ 
tal setup, similar to Nernst experiment, the Wiedemann-Franz 
law holds for both longitudinal and Hall conductivities (nor¬ 
mal and anomalous). For the parallel setup we hnd an ad¬ 
ditional B^ dependence of the Lorenz number arising from 
the chiral anomaly term v • In a previous theoretical 
worl^, both of these conclusions have been reported for a 
linearized WSM, and the violation of Wiedemann-Franz law 
in the parallel setup has been ascribed to the role of axion- 
electrodynamics because of the topological E • B term. We 
verify this violation of Wiedemann-Franz law in a lattice 
Hamiltonian, and it only depends on the presence of Berry 
curvature in a system, and therefore it is not an artifact of a 
linearized theory. 
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